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Abstract: 

The energy eigenvalues of the superintegrable chiral Potts model are determined by the 
zeros of special polynomials which define finite representations of Onsager's algebra. The 
polynomials determining the low-sector eigenvalues have been given by Baxter in 1988. In 
the Z3— case they satisfy 4-term recursion relations and so cannot form orthogonal sequences. 
However, we show that they are closely related to Jacobi polynomials and satisfy a special 
£> ■ "partial orthogonality" with respect to a Jacobi weight function. 
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£^ ■ 1 Introduction 



F.Y.Wu and Y.K.Wang |1[ were the first to consider the Potts model with chiral interaction 
terms. Their interest in this generalization arose from duality considerations, but the idea 
proved to be very fruitful in many respects: Ostlund Q and Huse ||] proposed the chiral 
Potts (CP) model for phenomenological applications: it allows to describe incommensurate 
phases using nearest neighbor interactions only. We give a few references [Q ||, H [?J from 
which the subsequent development can be traced, and turn directly to the superintegrable 
chiral 7Ln Potts quantum chain ||]. This is a particularly interesting model, because it 
provides some of the rare representations known for Onsager's algebra || and in this sense 
generalizes the Ising quantum chain (for N — 2 it is the Ising model). Integrability by 
Onsager's algebra entails that all eigenvalues of the hamiltonian are determined by the zeros 
of certain polynomials, which for the chiral Potts model were first derived by Baxter fl(i[| . 
Although the definition of Baxter's polynomials looks very simple, the properties of these 
polynomials turn out to be quite non-trivial and interesting Ji"l| . The main part of this note 
deals with the properties of these polynomials. They satisfy N + 1-term recursion relations, 
therefore for N > 2 they cannot form orthogonal sequences. However, as found recently 111]], 
several properties which characterize orthogonal polynomials are almost true for Baxter's 
polynomials (e.g. the zero separation property is true except for one extreme zero). 

We first recall the definitions of the superintegrable CP-hamiltonian and Onsager's algebra 
and then, following B.Davies fl2]| , we sketch how the formula for the energy eigenvalues 
emerges. We consider Baxter's polynomials and their recursion relations. Equivalent poly- 
nomials with their zeros in (—1, +1) for N = 3 are written in terms of a determinant. Their 
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expansion in terms of Jacobi polynomials gives the surprising result that many of the ex- 
pansion coefficients vanish, leading to the notion of "partial orthogonality". 

The hamiltonian defining the Z7v-superintegrable chiral Potts quantum chain ^] is: 

L N-l 



Here u> — e 27Tl ^ N and Zj and Xj are Z/v-spin operators acting in the vector spaces C N 
at the sites j — 1,2,... ,L (L is the chain length). The operators obey Zi Xj = 
XjZiUj SiJ ; Z^ = X? 1 = 1 and we assume Xl+i = X\ (periodic b.c). A convenient 
representation is [Xj)i >m = 5i, m +i mod N and (Zj)i, m — 8i. m uj m . For N > 3 the com- 
plex coefficients make the chain hamiltonian parity non- invariant. For N — 2 we get the 
Ising quantum chain. For fixed N there is only one parameter, the temperature variable 
k. Incommensurate phases arise due to ground state level crossings. TL commutes with the 
Z^r-charge Q = 1X^=1 ^i- We wr ite the eigenvalues of Q as ui® . Q — 0, 1, . . . , N — 1 labels 
the charge sectors of TL. 

We split TL into two operators writing TL = —h N(Ap + kAi). A remarkable property of TL 
is that Aq and A\ satisfy || the Dolan-Grady |T3] ] relations 

[Ao,[Ao,[A ,Ai]]] = \&[A a ,A{\- [A u [A u [A x , A ] ] ] = 16[A U A Q }, 

which are the conditions [0 for Ag and Ai to generate Onsager's algebra A, which is formed 
|)| from elements A m , G;, m G Z, I G N, I >m, satisfying 

[Ai,A m ]=4G l - m ; [G l ,A m ] = 2A m+l -2A m . l ; [G h G m ] = 0. (2) 

From (Q) , there is a set of commuting operators which includes TL: 

Q m = \{A m + A_ m + k{A m+1 +A- m+ x))\ [Qi, Q m ] = 0; Q = TL. 

To obtain finite dimensional representations of A we require the A m (and analogously the 
Gi) (l2[ [l5| to satisfy a finite difference equation: Sfe=-n ak Ak-i — 0. This is solved 
introducing the polynomial (the main object of the present paper): 

n 

T{z) = <*k z k+n (3) 

k— — n 

(from A the at are either even or odd in k). Now the A m and G m can be expressed in terms 
of the zeros Zj of T{z) and the set of operators , Hy. 

n n 

A m = 2 £ (zf E+ + zJ m Ej) ■ G m = Y. ( Z T ~ z J m ) H o ■ ( 4 ) 

From A these operators obey s/(2, C)-commutation rules: 

[E+E-] = 5 jk H k ; [H j7 E±\ = ±2 5 jk E±. 
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So A is isomorphic to a subalgebra of the loop algebra of a sum of sl(2, C) algebras. 
From the first of eqs.(^) we can express Ti in terms of the Zj and the operators Ej . Writing 
Ep = J x j ± iJy.j, then in a representation Z(n,s) characterized by the polynomial zeros 
Zi,...,z n and a spin-s representation of all the Jj, we get: 

n 

(Ao + kA,)^ = 2£{(2 + fc(^+z J - 1 ))4j+i(* i -*7 1 )43} 

n 

j'=i 

where J' X) j is a rotated 5'?7(2)-operator, and Cj = cosflj = + zj 1 ). 

For the CP-hamiltonians (Q) the spin representation turns out to be s — \ . Accordingly, all 
eigenvalues of ([!]) have the form 

E = -N (a + bk + 2 £™=i y/l + 2kcas9j + fc 2 ) , m 3 = ±±. (5) 
a and b are non-zero if the trace of Aq and Ai is non-zero. 



2 Baxter's polynomials 

No direct way is known to find the polynomials T{z) from the hamiltonian Ti. However, the 
invention of the two-dimensional integrable CP model 17 , which contains Ti as a special 
logarithmic derivative, and functional relations for its transfer matrix have enabled Baxter 
JTof to obtain the polynomials for the simplest sector of Ti, which at high-temperatures 
contains the ground state (the polynomials corresponding to all other sectors have been 
obtained subsequently in |l8, 19, 20[|). Here we shall consider only the simplest case. Baxter 



l0| finds that in terms of the variable t or s = t = (c — l)/(c + 1) (recall c = cos 8 of 
) these polynomials take the form 



P. 



."-1/1 



( s ) = iE T— 17 ) ^ t )~" Q - L : = ( Ar - !)( L + 0) mod ^ ( 6 ) 



Here Q denotes the Z^-charge sector. For Z 3 eq.(g) is, written more explicitly: 

P« ( s ) = { (t 2 + * + if + (t 2 + c 2 t + u) L + uj-Q (t 2 + ut + u; 2 ) L } . 

Due to their Zj\r-invariance t — ► ut, the depend only on s = t N . The degree of the 

Pq(s) in the variable s is bi,,Q — [((N — l)L — Q)/N] where [x] denotes the integer part 
of x. Considering sequences of these polynomials for fixed Q and L G N, we notice that 
the dimensions bL,Q do not always increase by one when increasing L by one: at every Nth 
step the dimension stays the same: e.g. the dimensions of the Pq—q for L mod N = and 
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L mod N = 1 coincide, see Table 1. 

The polynomials (g) have their zeros all on the negative real s-axis: TL is hermitian and so 
in (|^) we must have —1 < Cj < +1 which means negative Sj. We will prefer to deal mostly 
with equivalent polynomials in the variable c, defining 

n%\c) = (c + i) b ^pg\ s = §=i). (?) 

Our main concern in this paper is to learn about the properties of the H < q\c) or Pq j \s), 
e.g. whether these can be arranged into orthogonal sequences etc. We will find that the 
IIq^(c) are polynomials with quite remarkable properties. A number of special features of 
the IIq''' have been discussed recently pH . Here we give some more detailed results for the 

Z3-case. As the recursion relations for the Hq contain a lot of information, we now show 
how to obtain these. 



3 Recursion relations 



We start with the observation that the coefficients of the Pq(s) can be obtained from the 
expansion of (1 + t + t 2 + . . . + t N ~ v ) L , simply by taking every N— th term of the expansion, 
starting with the coefficient of t( Ar_1 ) L_< 3. More precisely, we claim that we can define the 



Pq 1 ^ by the decomposition 

{l + t + t 2 + ... + t N - l ) L = 



N 



l-t 
1 - t 



L N-l 
Q=0 



(8) 



demanding the P^ to depend on s — t N only. Proof: Insert (||) into (||) to get: 
l-t 1 



,N\L , N-l N-l ,.j(L+Q) 



l-t 



a-t N ) L - yy ± 

V ' N ^ ^ 1 



)=() .1=0 



LJO t) 



L ' 



Interchanging the Q- and j-summations we see that the Q-summation gives zero for j ^ 0, 

leaving only the j ' — term, which is (^) . 

Eq. (0) can now be used to obtain recursion relations: Write 



N-l 



N-l 



(l + t + t 2 + ... + t N - 1 ) taL ' Q p q\ s ) = t° L+± ' Q Pq +X \s) 

Q=0 Q=0 



(9) 



Comparing powers of t, e.g. for L mod N = this gives 



(P^ L+l) \ ( 1 



(L+l) 
(L+l) 



V p(L+l) I 
\ ^N-l I 



1 
1 

a 



1 \ 



V i 



(L) 
L 

(L) 



V P^ll J 



(10) 
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For L mod N = k replace Pq — > Pq_/c cyclically in both column vectors, keeping the 
same square matrix. Recursion relations not coupling P^ with different Q follow by the 
N— fold application of these relations, leading to N + 1-term recursion formulae. These can 
be transcribed into the corresponding formulae for the rig'''. For the Ising case Z2 these are 
of the Chebyshev type 

n£ +4 >-4 C n£ +2) + 4ii^ = o, (11) 

and so for N = 2 the ITq^ form orthogonal sequences. However, for Z3 we have JllJ (valid 
for all L > and all Q): 

n (L+9) _ 3 ( 9c 2 _ 5 ) n (i+6) + 48 n (£+3) _ 64 n (L) _ ( 12 ) 

These form 9 sequences, each labeled by (Lq, Q), where Q — 0, 1, 2 and L = 3 j + Lq 
where L = 0, 1, 2, j = 0, 1, 2, . . .. The degrees of the polynomials appearing in this relation 
increase by two from the right to the left, but since the recursion is four-term, not three- 
term, these are not orthogonal sequences [^l|[]. However, like ( p"l| ) also ( |l2"| ) are of the most 
simple typeR all coefficients are independent of L. 



4 Expansion in terms of Jacobi polynomials 

As the zeros of our are confined to and dense in the interval (—1, +1) we call this the 
basic interval like for orthogonal polynomials. Trying to determine (numerically) a weight 
function by the ansatz f\(l + c)"(l — cj^c'llg '(c) = for k < bQ t L fitting a and f3, we 
find (in the following we concentrate on the Z3-case) that there is an approximate solution 
very close to a = —(3 = i, but a and (3 come out to be slightly L and ^-dependent, in 
contrast to what is needed for orthogonality. However, for L — > 00 and small fixed k, the 
solutions converge towards a — —(3 = \- So the seem to be close to Jacobi polynomials 
(- 

P k 3 ' 3 , but can we formulate an exact relation valid for finite IS! Is there an exact property 
of the which replaces orthogonality? 

Numerical calculations (ll[] gave the surprising result that seemingly complicated can 
be written as a combination of just very few Jacobi polynomials, e.g. 

nf 2) = 3 n c 7 + 3 10 c 6 - 5 • 3 10 c 5 - 80919c 4 + 140697c 3 + 27459c 2 - 16839c - 1365 = 

TIE (If P7 3 ' 3 — ^5 3 3 ) ■ F° r polynomials ir(c) of degree n we useQ 



7T(C) = [W , TTi, 7T n ] = £ W k P^' ^(c) (13) 

fc=0 



1 If we consider Q = L mod 3, then we have only polynomials in c 2 , and the degrees 6_l,q are consecutive 
in powers of z = c 2 ("simple sets of polynomials") with integer coefficients. 

2 For higher TV we get similar N + 1-term relations, e.g. for Z4: 
n Q +16) - 4 ( 64c3 - 56c)l4 4) n^ +12) - 128(14c 2 - 17)n^ L+8) - 2048 C n^ +4) + 4096 n^ L) = 0. 

3 For Jacobi series as a generalization of Taylor series, see Ch.7 of Carlson ^3) . We use the standard 
normalization of Jacobi polynomials, see e.g. Rainville 122 1 
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Table 1: Examples of Z3— polynomials rig' (c) and their Jacobi-components defined in ([13]). 









Z 3 , Q = L + 1 


mod 3 




L 






n Q L, (c) 


2 -[2V3] n W( C ) 


3 


9c + 


3 








4 


27 c 2 


-18 c- 


-5 


[3 -21 £1 




5 


81 c 3 


+ 27c 2 


- 57c- 11 


[0 ^-1 




6 


3 5 c 3 


+ 81c 2 


- 135 c- 21 


[0,0,0,2f] 




7 


3 6 c 4 


- 2 • 3 5 


c 3 - 540c 2 + 270c + 43 


ro 27 171 729 
L°) 4 ' 14 ' 40 ' 


729] 
70 ' 


8 


3 7 c 5 


+ 3 6 c 4 


- 2754 c 3 - 702 c 2 + 711 c + 85 


[0,0,0,-f,0,^] 




9 


3 s c 5 


+ 3 7 c 4 


- 7290 c 3 - 1782 c 2 + 1593 c+ 171 


[0,0,0,-11,0,^] 




10 


3 9 c 6 


-2-3 8 


c 5 - 25515 c 4 + 14580 c 3 

+7965 c 2 - 3186 c -341 


ro 27 135 243 
L°) 4 ' 14 ' 20 ' 
3 7 
56' 


9477 
385 ' 

Jll 
308 J 



and define a scalar product with Baxter's variable t = ((1 — c)/(l + c)) 1 / 3 (see (||)) as the 
weight function (here we will not need to specify the normalization): 

(7rW|7r( 2 >> = [ +1 dc ( l^Y 7 ^ 1 ) (c)7r( 2 )(c) =- t dt , 6 * (c(s)) 7r^(c(s)). 

J-i Vi + cy J-oo (1-s) 2 

The second part of this definition shows that it makes sense also if we prefer to use polyno- 
mials in the variable s, and that it preserves the original Z3— symmetry. 

Since the IIq^ satisfy the recursion relations ( |l2[) they can be written as determinants of 
band matrices with a bottom line specifying the initial conditions (which are the lowest L 
polynomials). Omitting the bottom line, we define j x j— band matrices and their determi- 
nants Rj = det j band( [64, 48, 3p, 1, 0], j)\ , e.g. 



i?6 — 



3p 


1 














48 


3p 


1 











64 


48 


3p 


1 











64 


48 


3p 


1 











64 


48 


3p 


1 











64 


48 


3 P 



where p = 9c 2 — 5, so that the polynomials Rj depend only on c 2 . Now we get the nine 
sequences of the IIq^ as linear combinations of the Rj which satisfy appropriate initial con- 
ditions. Abbreviating Qj = Rj + 8i?,_i these are found to be: 
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n o 3j) = § (Qj ~ 8Qj-i + 16Q 3 --a) ; n^' 5 = 9c%_i ± 3Qj_i; 

2 

n (3j+i) = ± 18cRj l + q. + 2 q . nf i+1) = Q i - 4Q J _ 1 ; (14) 

2 

n ( 3j+2) = 3c(Rj - 4R,-_i) ± (Qj - 4Q,_ 1 ); n^' +2) = 3Q,. 

i 

For j = 1, 2 use i?o = 1; R-i = R-2 = 0. From ( |l4| ) we see that only two of the 9 sequences 
are independent. There are relations like e.g. IIg ' +1 — = 2IIj . 

To get the Jacobi-expansion of the IIq^ , we only need to expand Rj and cRj . Using Jacobi- 
components defined in (|l3|), from explicit calculation (for j < 36), we find that for k < j 
(only there) the j-dependence of the (Rj) k obeys the simple rule: 

(R,) k = (-3) fe k\ (2k + 1) {(-8Y a k + 4> r k ) ; 

1 _ 2(-3) m _ _ 2(-3) m ~ 2 

3nn=o(3«+l) 3n'= m (2« + l) Ul= m (^) 

It follows that the a k do not contribute to the k < j-components of Qj, and, using the 

(A _i) 

recursion relations for the 3 ' 3 (c) , we conclude that for fc < j we have 
i(Q,) fc = -(cR j ) k = \(c 2 R ]+1 ) k = 4^(-3) fc kl (2k + l)r k . 
It further follows that 

(u^) k = (n[^) k ^(u^) k = (u^) = (irT +2 >) =o for k<j. 

i k i k 

All fc < j components of II2 and H^ +2 ^ are proportional to (Qj) k - 

We get zero overlap between a polynomial ITq^ of degree &l,q with all polynomials Hq', ^ 
which have at least &l,q vanishing low-fc components (these can only be polynomials which 
have bL'jQ 1 > 2&l,q)- One of many such relations is e.g. 

(ng j) I ng/' } ) = for Q = Q'= 1 and 2j < f - 1. 
This property may be called "partial orthogonality". 

Further rules, this time valid for all fc, regard the vanishing of many components for partic- 
ular linear combinations: Defining 

_r~(+) „.(+) „.(+) 1. 



2j-l • 



Q+ = = 9c^_i +3Q,_! - [a +; , a™, . .. , a 

we have checked up to j = 30 that a k — for k < j and that all even (odd) k Jacobi 
components of Qr) (Q^) vanish. So we conjecture for all j, j' 

(Q ( ?\Q { P) = 0. 

One can check some special cases of these results in Table 1. It has been found numerically 
[ p"l| that a similar partial orthogonality appears also for the Z^r-Baxter polynomials with 
N = 4, 5, 6. For even values of N further relations emerge, but these will not be discussed 
here. 
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5 Conclusion 



The polynomials which play a central role for the calculation of the energy eigenvalues of 
the superintegrable Z3— chiral Potts model are found to be related to Jacobi polynomials in 
a very peculiar way. Many integrals giving the Jacobi-coefficients of Baxter's polynomials 
are found to vanish. These observations should have a deeper group-theoretical background, 
but the underlying symmetry is not yet clear to us. By reducing the formulation of the 
problem to some basic facts, the present analysis tries to prepare the ground for clarifying 
the symmetry involved. 
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